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Abstract 

Discovery potentials for extra neutral interactions at the Large Hadron Collider 
in forthcoming experiments are analyzed using resonant leptoproduction. For this 
purpose we use high precision next-to-next-to-leading order (NNLO) determinations 
of the QCD background in this channel, at the tail of the Drell-Yan distributions, 
in the invariant mass region around 0.8 < Q < 2.5 TeV. We focus our analy- 
sis primarily on a novel string-inspired Z', obtained in left-right symmetric free 
fermionic heterotic string models and whose existence at low energies is motivated 
by its role in suppressing proton decay mediation. We analyze the parametric de- 
pendence of the predictions and perform comparison with other models based on 
bottom up approaches, that are constructed by requiring anomaly cancellation and 
enlarged Higgs structure. We show that the results are not particularly sensitive to 
the specific charge assignments. This may render quite difficult the extraction of 
significant information from the forward-backward asymmetries on the resonance, 
assuming that these are possible due to a sizeable width. The challenge to discover 
extra (non anomalous) Z' in this kinematic region remains strongly dependent on 
the size of the new gauge coupling. Weakly coupled extra Z' will not be easy to 
identify even with a very good theoretical determination of the QCD background 
through NNLO. 



1 



1 Introduction 



The search for neutral currents mediated by extra gauge bosons (Z') at the Large Hadron 
Collider will gather considerable attention in the next few years [1]. Additional Abelian 
gauge interactions arise frequently in many extensions of the Standard Model, like in 
left-right symmetric models, in Grand Unified Theories (GUTs) and in string inspired 
constructions [lj. It has also been suggested that the existence of a low scale Z' may 
account for the suppression of proton decay mediating operators in supersymmetric theo- 
ries and otherwise [21 El S] • Abelian gauge structures may also play a considerable role in 
fixing the structure of the flavor sector, for instance in pinning down the neutrinos mass 
matrix. Anomaly cancellation conditions, when supported also by an extended Higgs 
and fermion family structure - for instance by the inclusion of right-handed neutrinos - 
may allow non-sequential solutions (i.e. charge assignments which are not proportional to 
the hypercharge) that are phenomenologically interesting and could be studied by ATLAS 
and CMS. Furthermore, within left-right symmetric models, and their underlying SO (10) 
embedding, the global baryon minus lepton number (B — L) of the Standard Model is 
promoted to a local symmetry. Abelian gauge extensions are therefore among the most 
well motivated extensions of the Standard Model. For these reasons, the identification 
of the origin of the extra neutral interaction in future collider experiments will be an 
important and challenging task. In particular, measurements of the charge asymmetries 
- both for the rapidity distributions and for the related total cross section - and of the 
forward-backward asymmetries, may be a way to gather information about the structure 
of these new neutral currents interactions, although in the models that we have studied 
this looks pretty difficult, given the low statistics. 

As an extra Z' is common in model building, the differences among the various con- 
structions may remain unresolved, unless additional physical requirements are imposed 
on these models in order to strengthen the possibility for their unique identification. In 
this work we analyze the potential for the discovery of an extra Z' arising in a specific 
string construction, which is motivated not only by an anomaly-free structure, as in most 
of the bottom-up models considered in the previous literature, but with some additional 
requirements coming from an adequate suppression of proton decay mediation. Bottom 
up approaches based only on anomaly cancellation are, in this respect, less constraining 
compared to models derived either from a string construction or from theories of grand 
unification (GUTs) and can only provide a basic framework within which to direct the 
experimental searches. At the same time the search for extra neutral interactions has to 
proceed in some generality and be unbiased, looking for resonances in several complemen- 
tary channels. In this work we will investigate the relation between more constrained and 
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less constrained searches of extra neutral gauge bosons by choosing as a channel lepto- 
production and proceed with a comparison of some proposals that have been presented 
in the recent literature. Our main interest is focused around an extra Z' which has been 
derived using the free fermionic formulation of string theory in a specific class of left -right 
symmetric string models. The new abelian structure is determined not just as an attempt 
to satisfy some additional physical requirements, on which we elaborate below, but is 
naturally derived from a class of string models which have been extensively studied in 
detail in the past two decades [HI El E|. 

Our paper is organized as follows: in section [2] we discuss the origin of Z' in heterotic- 
string models. We discuss in some details the origin of the charge assignment under the 
Z' , which is motivated from proton decay considerations and differs from those that have 
traditionally been discussed in the literature. Then we move to define the conventions in 
regard to the charge assignments and the Higgs structure of the models that we consider, 
which are characterized by a gauge structure which enlarges the gauge group of the Stan- 
dard Model by one extra U(l). Our numerical analysis of the invariant mass distributions 
for leptoproduction is performed by varying both the coupling of the extra U(l) and the 
mass of the new gauge boson. The dependence on these parameters of the models that 
we discuss are studied rather carefully in a kinematic region which can be accessed at the 
LHC. We compare these results with those obtained for a group of 4 different models, 
introduced in [TO], for which we perform a similar analysis using leptoproduction. From 
this analysis it is quite evident that the search for extra neutral currents at the LHC 
is a rather difficult enterprise in leptoproduction, unless the coupling of the new gauge 
interaction is quite sizeable. 



2 Heterotic— string inspired Z' 

Phenomenological string models can be built in the heterotic-string or, using brane con- 
structions, in the type I string. The advantage of the former is that it produces states in 
spinorial representations of the gauge group, and hence allows for the SO (10) embedding 
of the matter spectrum. The ten dimensional supersymmetric heterotic-string vacua give 
rise to effective field theories that descend from the E 8 x E 8 or 5*0(32) gauge groups. The 
first case gives rise to additional Z's that arise in the SO (10) and Eq extensions of the Stan- 
dard Model, and are the cases mostly studied in the literature pQ. A basis for the extra Z' 
arising in these models is formed by the two groups U(l) x and U(l)^ via the decomposition 
E 6 -> 50(10) x £7(1)^ and 50(10) -> SU(5) x U(l) x pQ. Additional, flavor non-universal 
t/(l)'s, may arise in heterotic Eg x E 8 string models from the U(l) currents in the Cartan 
subalgebra of the four dimensional gauge group, that are external to Eq. Non-universal 
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Z's typically must be beyond the LHC reach, to avoid conflict with Flavor Changing Neu- 
tral Currents (FCNC) constraints. Recently [I] a novel Z' in quasi-realistic string models 
that do not descend from the heterotic E$ x E 8 string has been identified. Under the new 
£7(1) symmetry left-handed components and right-handed components in the 16 spinorial 
5*0(10) representation, of each Standard Model generation, have charge —1/2 and +1/2, 
respectively. As a result, the extra U(l) is family universal and anomaly free. It arises in 
left-right symmetric string models |9|, in which the 50(10) symmetry is broken directly 
at the string level to SU(3) x U(1) B -l x SU(2) l x SU{2) r x U(1) z > x U(l) n x hidden 
[9]. The U(l) n are flavor dependent U(l)s that are broken near the string scale. The 
Standard Model matter states are neutral under the hidden sector gauge group, which in 
these string models is typically a rank eight group. It is important to note that the fact 
that the spectrum is derived from a string vacuum that satisfies the modular invariance 
constraints, establishes that the model is free from gauge and gravitational anomalies. 
The pattern of U(l)z' charges in the quasi-realistic string models of ref. [9] does not arise 
in related string models in which the 50(10) symmetry is broken to the SU(5) x U(l) 
[§], the 50(6) x 50(4) 0, or SU(3) x SU(2) x U(l) 2 [H], subgroups. The reason for the 
distinction of the left-right symmetric string models is the boundary condition assignment 
to the world-sheet free fermions that generate the 50(10) symmetry in the basis vectors 
that break the 50(10) symmetry to one of its subgroups. The world-sheet fermions that 
generate the rank eight observable gauge group in the free fermionic models are denoted 
by {^r''"' 5 , 77 1 ' 2 ' 3 }, where if) 1 ''"' 5 generate an 50(10) symmetry, and fj 1,2,3 produce three 
U(l) currents^. Additional observable gauged U{1) currents may arise at enhanced sym- 
metry points of the compactified six dimensional lattice. The 50(10) gauge group is 
broken to one of its subgroups SU(5) x U(l), 50(6) x 50(4) or SU(3) x SU{2) x U(l) 2 

by the assignment of boundary conditions to the set ip\" 5 : 

2 

1. &{<P"V' 2 ' 3 } = =► su ^) >< x u(i)\ (1) 

2. ■ 5 r7 1 ' 2 ' 3 } = {11100000} 50(6) x 50(4) x U(l) 3 . 

To break the 50(10) symmetry to§ SU(3) C x SU(2) L x U(l) c x U(1) L both steps, 1 
and 2, are used, in two separate basis vectors. The breaking pattern 50(10) — > SU(3)c x 
SU (2) l x SU {2)r xU(1) b-l is achieved by the following assignment in two separate basis 
vectors 

1. V' 2 ' 3 } = {11100000} 50(6) x 50(4) x U{1) 3 , (2) 

2. K^-V' 2 ' 3 } = ^ SU $)c x U{l)c x SU(2) L x SU(2) R x U(l) 3 



^^for reviews and the notation used in free fermionic string models see e.g. [S] and references therein. 
2 U(1) C = \U(X) B -l\ U(1) l = 2U{l) T . 
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The distinction between the symmetry breaking patterns in eq. (CD) and eq. (j2J) is 
with respect to the charges of the Standard Model states under the three flavor depen- 
dent U(l) symmetries (1)1,2,3 that arise from the three world-sheet fermions 77 1 ' 2 ' 3 . In 
the free fermionic models, the states of each Standard Model generation fit into the 16 
representation of 5*0(10), and are charged with respect to one of the three flavor U(l) 
symmetries. For the symmetry breaking pattern given in eq. (JTJ) the charge is always 
+1/2, i.e. 

Q 3 (16 = {Q,L,U,D,E,N}) = +^ (3) 
whereas for the symmetry breaking pattern in eq. (j2J) the charges are 



Q](Ql, Ll) = -- 
Q j (Q R ={U,D},L R = {E,N}) = + l - (4) 

As a result in the models admitting the symmetry breaking pattern eq. (pP) the combina- 
tion 

I7(l) f = tf(l)i + tf(l) 2 + tf(l) 3 . (5) 

is anomalous, whereas in the models admitting the symmetry breaking pattern (TSJ) it is 
anomaly free. The distinction between the two boundary condition assignments given in 
eqs. ([T]) and ([2]), and the consequent symmetry breaking patterns, is important for the 
following reason. Whereas the first is obtained from an N = 4 vacuum with E$ x E$ 
or 5*0(16) x 50(16) gauge symmetry, arising from the {■0 1 ''"' 5 , fj 1,2 ' 3 ^ 1 ''"' 8 } world-sheet 
fermions, which generate the observable and hidden sectors gauge symmetries, the second 
cannot be obtained from these N = 4 vacua, but rather from an N = 4 vacuum with 
50(16) xE^xEj gauge symmetry, where we have included here also the symmetry arising 
from the compactified lattice at the enhanced symmetry point. The important fact from 
the point of view of the Z' phenomenology in which we are interested is that the first case 
gives rise to the type of string inspired Z' that arises in models with an underlying E e 
symmetry. Whereas the Eq may be broken at the string level, rather than in the effective 
low energy field theory, the crucial point is that the charge assignment of the Standard 
Model states is fixed by the underlying Eq symmetry. The entire literature on string 
inspired Z' studies this type of Eq inspired Z' . The second class, however, is novel and 
has not been studied in the literature. In this respect it would be interesting to examine 
how the symmetry breaking pattern ([2D and the corresponding charge assignments (111) 
can be obtained in heterotic orbifold models in which one starts from a ten dimensional 
theory and compactifies to four dimensions, rather than starting directly with a theory 
in four dimensions, as is done in the free fermionic models. This understanding may 
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highlight the relevance of ten dimensional backgrounds that have thus far been ignored 
in the literature. From the point of view of the Z' phenomenology, which is our interest 
here, the crucial point will be to resolve between the different Z' models and the fermion 
charges, which will reveal the relevance of a particular symmetry breaking pattern. 

The existence of the extra Z' at low energies, within reach of the LHC, is motivated by 
proton longevity, and the suppression of the proton decay mediating operators [21 El H]- 
The important property of this Z 1 is that it forbids dimension four, five and six proton 
decay mediating operators. The extra U(l) is anomaly free and family universal. It allows 
the fermions Yukawa couplings to the Higgs field and the generation of small neutrino 
masses via a seesaw mechanism. String models contain several U(l) symmetries that 
suppress the proton decay mediating operators [3]. However, these are typically non- 
family universal. They constrain the fermion mass terms and hence must be broken at 
a high scale. Thus, the existence of a U(l) symmetry that can remain unbroken down 
to low energies is highly nontrivial. The £7(1) symmetry in ref. [HI H] satisfies all of 
these requirements. Furthermore, as the generation of small neutrino masses in the string 
models arises from the breaking of the B — L current, the extra U(l) allows lepton number 
violating terms, but forbids the baryon number violating terms. Hence, it predicts that 
i?-parity is violated and its phenomenological implications for SUSY collider searches 
differ substantially from models in which i?-parity is preserved. The charges of the 
Standard Model states under the Z' are displayed in table [91 Also displayed in the table 
are the charges under U(l)^ = Uc — Ul, which is the Abelian combination of the Cartan 
generators of the underlying SO (10) symmery that is orthogonal to the weak hypercharge 
U(l)y- The charges under the U(l) combination given in eq. ((Sj) are displayed in table 
[9] as well. These two Z7(l)'s are broken by the VEV that induces the seesaw mechanism, 
and the combination 

U(l) z , = \u(l\, - U(l\ (6) 

is left unbroken down to low energies in order to suppress the proton decay mediating 
operators. The charges of the Standard Model states under this U(l)z> are displayed in 
table El 

3 The interactions for U(\)z' 

In this section we fix our conventions and describe the structure of the new neutral sector 
that we are going to analyze numerically in leptoproduction afterwards. The notations 
are the same both in the case of the string model and for the other models that we will 
investigate. We show in © the field content of the string model obtained within the free 
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fermionic construction discussed above. Of the 3 extra £7(1), we will decouple the two 
gauge bosons denoted by (, (' and keep only the Z'. The assumption of decoupling of 
these extra components are realistic if they are massive enough (> 5 TeV) so to neglect 
their influence on the lowest new resonance. We have chosen a mass Mz> around 0.8 TeV. 
We recall that a reasonable region where the new extra gauge boson have a chance of 
being detected is below the 5 TeV range. 

The fermion-fermion-Z' interaction is given by 

^ZfgJ-ffZl, (7) 
/ 

where / = e R , l 3 L , u R , d R , q 3 L and q 3 L = (u 3 L ,d 3 L ) ,l 3 L = {v 3 L) e 3 L ). The coefficients z u ,z d are 
the charges of the right-handed up and down quarks, respectively, while the z q coefficients 
are the charges of the left-handed quarks. g z is the Z' coupling constant. We can write 
the Lagrangean for the Z'-lepton-quark interactions as follows 



(8) 



with j being the generation index. The low energy spectrum of the model, as discussed 
above, is assumed to be the same for the other models that we analyze in parallel. As 
shown in ([9]) the field content of the model is effectively that of the Standard Model 
plus 1 additional Higgs doublet. The extra scalars 0, and (h and the right handed 
components iV# and iV# are assumed to decouple. In this simplified framework, the 
structure of the vertex 
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(9) 



is the following 
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where V M denotes generically the vector boson. In the Standard Model (SM) 
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(10) 
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(11) 



We need to generalize this formalism to the case of the Z'. 

Our starting point is the covariant derivative in a basis where the three electrically- 
neutral gauge bosons B Y , B% are 

■ 9y 



D, 



d, - ig (W^T 1 + W 2 T 2 + WjT 3 ) - ifYB Y - i^B* 



(12) 



and we denote with g,gy,g z the couplings of SU(2), U(l)y and U(l) z , with tan6V = 
gyjg. After the diagonalization of the mass matrix we have 

/ A u \ ( sin0 w cos0 w \ ( Wl\ 



\ Z 'J 



cos 6 W — sin 9 W e 
y — e sin £ sin @w 1 



V B l ) 



(13) 
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where e is defined as a perturbative parameter 



Then we define 



_ 5M 2 ZZ , 
6 ~ Mf, - Mf 

Mf = <r 



4 cos 2 



4+^ 2 )[l + 0(e 2 )] 



m| = f (4A + 4 2 < + 44) [i + o(^ 2 )] 



5M 2 ZZ , = 



99z 



4 cos 9w 



I 2 2 , 2 2 \ 



(14) 



9 = — 



9y = 



(15) 



sin 9w cos 6^ 

and we construct the W ± charge eigenstates and the corresponding generators T* 1 as 
usual 

W 1 =f iW 2 



± 



W - -= 

, _ Ti ± iT 2 
V2 ' 



(16) 



with the rotation matrix 




/ smd w (1+e 2 ) cos gyy ^_ cos 6> w 
' 1+e 2 1+e 2 fc 1+e 2 

cos 9 w (l+e 2 ) 

1+e 2 ~ 1+e 2 

e 1 



sin gyy sin_9w_ 

t 1+£ 2 







1+e 2 



1+e 2 



(17) 



from the interaction to the mass eigenstates. Substituting these expression in the covariant 
derivative we obtain 



D, 



Y 



- gT 3 sin 9 W + g Y cos 9 w -)-ig (W~T~ + W+T+) 



■iZ,, \ g cos 9 W T ?> - g Y sin 8 W ^ + g z e^ 



Y 



- iZ 'm [ ~9 cos 9 w T 3 e + g Y sin 9 W ^—e + g z | 



(18) 



where we have neglected all the 0(e 2 ) terms. Sending g z — > and e — > we obtain the 
SM expression for the covariant derivative. The next step is to separate left and right 
contributions in the interactions between the fermions and the Z' boson. Hence for the 
quarks and the leptons we can write an interaction Lagrangean of the type 

C int = Q J L NfYQ J L z, + q{n z l ' yoiz; + u R N^u R Z, 
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Figure 1: Plot of the LO, NLO and NNLO cross section for the free fermionic model with 
M z > = 800 GeV. 



+d j R N!^d j R Z, + u^N^Yu^ + d^N^d^ 
+Q 3 L Nl^Q{A, + u^Nl^u^A, + d^Nj^d^A, 
+V L NIYV L A, + e^Nl^A, 
+V l N z LMj1 %Z, + V L N%jrfff L ^ 



where for the quarks we have 



N. 



( Y L 

l 3 = -M9 cos 9 w Tz - g Y sin W ~^ 

I' 3 = -il-g cos 9 w T^e + g Y sin 9 w —e + g z — 



N Z d , R = -i 



Y u >R 2 U >R 
-g Y sin + 9z£— 



z d,R 



-g Y sin 9w—^~ + 9 z e— 



(19) 



(20) 



and similar expressions for the leptons. We rewrite the vector and the axial coupling of 
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the Z and Z' bosons to the quarks as 

-ig 1 



z?V^ = — - 



9 



r 



T9v = 



9 



4c 



+ *4( y + f ) + f ^ + ) 



4c 



c w 2 



(21) 



where j is an index which represents the quark or the lepton and we have set sin#yK 
s w , cos^vk — c w for brevity. 

The decay rates into leptons for the Z and the Z' are universal and are given by 



Y(Z -> 11) 

T{Z -> Vi^i) 
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[ 1 + ^ + 1 . 409 ^- 12.77^ 

7T 7T 2 

where i = u,d,c, s and Z = Z, Z' . 

For the Z' and Z decays into heavy quarks we obtain 

N r a, 



(22) 
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(23) 



The total hadronic widths are defined by 

V z = r(Z — > hadrons) = r(Z — > V'iV'i) 
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Figure 2: Plot of the LO, NLO and NNLO cross section for the free fermionic model with 
M z > = 800 GeV in the TeVs region. 



Y z 1 = T(Z' — > hadrons) = S ^^T{Z' — > ipiipi 



(24) 



where we refer to hadrons not containing bottom and top quarks (i.e. i = u,d,c,s). 
We also ignore electroweak corrections and all fermion masses with the exception of the 
top-quark mass, while we have included the relevant QCD corrections. Similarly to pH] 
we have considered only tree level decays into fermions, assuming that the decays into 
particles other than the SM fermions are either invisible or are negligible in their branching 
ratios, then the total decay rate for the Z and Z' is given by 

T z = Yl T(Z^^ l ) + T(Z^bb) + 'ZT(Z^ll) + ?,T(Z^v l v l ) 

i=u,d,c,s 

T z <= F(Z' ^ ^ifji) + T(Z' ^bb) + T(Z' ^tt) + 3T(Z' ^U) + 3T(Z' ^ m)- 

(25) 



i=u,d,c,s 



We also recall that the point-like cross sections for the photon, the SM Z and the 
new Z' gauge boson are written as 

Ann 2 1 



3Q 4 iV c 



a z (Q 2 ,M 2 z ) 



z^u 



4Mz sin 2 9 W cos 2 6 W N C (Q 2 - M 2 ) 2 + M 2 T 2 
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Figure 3: Free fermionic model at the LHC, tan/3 = 40 
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(26) 



where Nc is the number of colours, and 



4M Z , sin 2 9 W cos 2 6 W N C (Q 2 - M 2 ,) 2 + M 2 T|, 
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(Q 2 - M|)(Q 2 - Mf,) + M z r z M z ,r z , 



96 sin 4 W cos 4 w iV c [(Q 2 - M 2 ,) 2 + Mf,r|,] \(Q 2 - Mf ) 2 + M 2 Y\ 



(27) 



The contributions such as Z, 7 and similar denote the interference terms. At LO (or 
leading order) the process proceeds through the qq annihilation channel and is 0(1) in 
the strong coupling constant a s . The NLO (or next-to-leading order) corrections involve 
virtual corrections with one gluon exchanged in the initial state and real emissions in- 
volving a single gluon, which is integrated over phase space. These corrections are 0(a s ) 
in the strong coupling. The change induced by moving from LO to NLO amounts to ap- 
proximately a 20 to 30 % in the numerical value of the cross section that we consider. At 
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the highest accuracy, we use in our analysis partonic contributions with hard scattering 
computed at NNLO, or 0(a 2 ). At this order typical real emissions involve 2 partons in the 
final state - which are integrated over their phase space- and two-loop virtual corrections 
at the same perturbative order. The cross section for the invariant mass distributions 
factorizes at a perturbative level in terms of a NNLO (next-to-next-to- leading, or 0(a 2 )) 
contribution Wy (which takes into account all the initial state emissions of real gluons and 
all the virtual corrections) and a point-like cross section. The computation of Wy can be 
found in [llj to which we refer for more details. A similar factorization holds also for the 
total cross section if we use the narrow width approximation. At NLO (next-to-leading 
order, or 0(a s )). The colour-averaged inclusive differential cross section for the reaction 
p + p — > l\ + l 2 + A, is given by 

<h ra v (Q 2 ,M 2 v )Wy(r,Q 2 ) r = (28) 



dQ2 .-V^ ,~ V/ ..y K .,^ , s 

where all the hadronic initial state information is contained in the hadronic structure 
function which is defined as 

W v (t,Q 2 ) = V" / dxi I dx 2 / dx5(r - xx 1 x 2 )PDY j (x 1 ,x 2 , ^F)Aij(x,Q 2 , /jp) , 
itj Jo Jo Jo 

(29) 

where the quantity PDYAxi, x 2 ,fi 2 F ) contains all the information about the parton distri- 
bution functions and their evolution up to the \ip scale, while the functions Ay (a;, Q 2 , fi 2 F ) 
are the hard scatterings. This factorization formula is universal for invariant mass dis- 
tributions mediated by s-channel exchanges of neutral or charged currents. The hard 
scatterings can be expanded in a series in terms of the running coupling constant a s (fi 2 R ) 
as 

oo 

Ay (*, Q\ /4) = <(^) A £- Q 2 ' PF, A) ■ (30) 

n=0 

In principle, factorization and renormalization scales should be kept separate in order 
to determine the overall scale dependence of the results. However, as we are going to 
show, the high-end of the Drell-Yan distribution is not so sensitive to these higher order 
corrections, at least for the models that we have studied. 



4 Numerical Results 

In our analysis we have decided to compare our results with a series of models introduced 
in [TO] - We refer to this work for more details concerning their general origin. We just 
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mention that the construction of models with extra Z' using a bottom-up approach is, in 
general, rather straightforward, being based mostly on the principle of cancellation of the 
gauge cubic U{1)\, and mixed anomalies. One of the most economical ways to proceed is 
to introduce just one additional SU(2)w Higgs doublet and an extra scalar (weak) singlet, 
as in P2], and one right-handed neutrino per generation in order to generate reasonable 
operators for their Majorana and Dirac masses. However, more general solutions of the 
anomaly equations are possible by enlarging the fermion spectrum and/or enlarging the 
scalar sector [13]. In [TU] the scalar sector is enlarged with 2 Higgs doublets and one 
(weak) scalar singlet. 

Anomalous constructions, instead, require a different approach and several phenomeno- 
logical analysis have been presented recently [TU [T5J, [T7] that try to identify the sig- 
nature of these peculiar realizations. In the anomalous models, due to the absence of the 
non-resonant behaviour of the s-channel (at least in the double prompt photon produc- 
tion), the chiral anomaly induces a unitarity growth which should be present in correlated 
studies of other channels [T7j. For non anomalous Z' the phenomenological predictions 
are, as we are going to show, rather similar for all the models - at least in the mass invari- 
ant distributions in Drell-Yan - and the possibility to identify the underlying interaction 
requires a careful study of the forward-backward and/or charge asymmetries [18] . This 
is not going to be an easy task at the LHC, given the size of the cross section at the tail 
of the invariant mass distribution, the rather narrow widths, and given the presence of 
both theoretical and experimental errors in the parton distributions (pdf's), unless the 
gauge coupling is quite sizeable (0(1)). We refer to [19] for an accurate analysis of the 
experimental errors on the pdf's in the case of the Z peak. It has been shown that the 
errors on the pdf's are comparable with the overall reduction of the cross section as we 
move from the NLO to the NNLO. 

These source of ambiguities, known as experimental errors, unfortunately do not take 
into consideration the theoretical errors due to the implementation of the solution of the 
DGLAP in the evolution codes, which amount to a theoretical uncertainty [20]. Once 
all these sources of indeterminations are combined together, the expected error on the 
Z peak is likey to be much larger than 3 %. Given the large amount of data that will 
accumulate in the first runs (for Q = Mz), which will soon reduce the statistical errors on 
the measurements far below the 0.1 % value, there will be severe issues to be addressed 
also from the theoretical side in order to match this far larger experimental accuracy. The 
possibility to use determinations of the pdf's on the Z peak for further studies of the Z' 
resonances at larger invariant mass values of the lepton pair, have to face several additional 
issues, such as the presence of an additional scale, which is Q = Mz', new respect to the 
Q = Mz scale used as a benchmark for partonometry in the first accelerator runs. We 
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Figure 4: Free fermionic model at the LHC, g z — 0.1 

remind that logarithms of these two scales may also play a role especially if Mz> is far 
larger than Mz- With these words of caution in mind we proceed with our exploration 
of the class of models that we have selected, starting from the string model and then 
analizing the bottom- up models mentioned above [TU]. These are studied in the limit 
= zh 2 = 0, with the mass of the extra Z' generated only by the extra singlet scalar 0. 
In the string model, as one can see from ([9]), only the two Higgses Hjj and Hp contribute 
to the mass of the new gauge boson. The differences between these two types of models 
are, however, not relevant for this analysis, since the mass of the extra gauge boson is 
essentially a free parameter in both cases. 

The set of pdf's that we have used for our analysis is MRST2001 [2JJ, which is given 
in parametric form, evolved with CANDIA (see [22] ). The models analyzed numerically 
are the free fermionic one, "F", discussed in the previous sections, and the U B — L", "q 
+u" , "10 + 5" and "d-u", using the notations of [TO] . 

Our results are organized in a series of plots on the various resonances and in some 
tables which are useful in order to pin down the actual numerical value of the various 
cross sections at a given invariant mass. 
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4.1 M z > = 0.8 TeV 



We show in Fig. [T]a plot of the Z' resonance around a typical value of 800 GeV for the FF 
model and the SM. The coupling of the extra neutral gauge boson is taken to be 0.05, with 
tan/3 = 10. We remark that the dependence of the resonance on this second parameter is 
negligible. In fact the relevant parameters are the coupling constant gz and the mass Mz>- 
Notice that the width is very narrow (~ 1 GeV) and basically invisible in an experimental 
analysis. Neverthless it is, at least theoretically, useful to try to characterize the signal 
and the background even in this (and other similar) not favourable cases. 

Assuming an integrated luminosity of 100 fb^ 1 /y after the first 3 years at the LHC 
(per experiment), we would expect 10 background events versus a signal of approximately 
30 events. Notice that LO, NLO and NNLO determinations are, essentially, coincident 
for all the practical purposes. 

In Fig. [2] we show the tail of the distribution for a run with Mz> = 800 GeV, where 
we have just modified tan/3 and we have increased the coupling to gz = 0.1. For Q 
around 1.2 TeV the determinations of the cross section in the FF and SM models are 
basically overlapping as we move from LO to NLO and NNLO. The LO determination 
in the SM moves up toward the FF result as we increase the perturbative order. Also in 
this case, given the small size of the cross section (^ 10~ 2 fb) the possibility to resolve 
these differences experimentally is remote. In Fig. [3] we vary the coupling constants of 
the extra U(l) from a very small value gz = 0.05 up to gz = 0.2. The only variation in 
the result is due to the width that increases from 1 to approximately 3-4 GeV's. Here we 
have chosen tan /3 = 40, and, as shown in Fig. H] there is essentially no variation on the 
shape of the resonance due to this variable. In Figs. and [U] we perform a comparative 
study of all the models and the SM background for a resonance mass of 800 GeV. There 
are only minor differences between the 4 bottom-up models and the FF model. The FF 
model shows a resonance curve which sits in the middle of all the determinations but is, 
for the rest, overlapping with the other curves. The " B — L" model, in all the cases, 
shows a wider width among all, with the "g + u n model quite similar to it. The u d — u n 
model has the narrowest width. This feature is particularly obvious from Fig. [7] where 
the result is numerically smoothed out by the increased value of the coupling, which is 
now doubled compared to Fig. El 

4.2 M z > = 1.2 TeV 

We illustrate in the next 3 figures our results for the various models for Mz> = 1.2 TeV. 
Fig. [H] shows the behaviour of the cross section for this new mass value with gz = 0.1, 
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Figure 5: Free fermionic model at the LHC, tan/? = 40 and g z = 0.05 

and the corresponding result for the SM case. The QCD corrections are very small and 
it is likely that the only role of these corrections, at these large Q values, is to stabilize 
the dependence of the perturbative series from the factorization/renormalization scales. 
In our case we have chosen, for simplicity fip = Hr = Q, where fiR and fip are the 
renormalization and factorization scale, respectively. The separation of this dependence 
can be done as in [20], by relating the coupling constants at the two scales (/ip, /Jr). 

This separation, in general, needs to be done both in the hard scattering and in the 
evolution. A zoom of the resonance region is shown in Fig. [HI which shows that the 
reduction of the signal is by a factor of 10 compared to the case of Mz> = 0.8 TeV. This 
drastic reduction of the cross section is one of the reason why the search of extra neutral 
currents, if these are mediated by new gauge bosons of mass above the 1 TeV range, may 
take several years of LHC luminosity to be performed, unless the new gauge coupling 
is larger. As shown in Fig. [TDl as we move away from the resonance region, the SM 
background and the FF result overlap. An interesting feature is that the K-factors for 
the SM result are much larger than for the FF case, especially as we move from LO to 
NLO. At NNLO both curves, however, overlap. 

We show, in Fig. [TT] a plot of the shape of the resonance region for Mz> = 2.5 TeV. 
The width is very narrow (2 GeV) and the size of the cross section down by a factor of 100 
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Figure 6: Free fermionic model at the LHC, tan/3 = 40 and g z — 0.1 
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Figure 7: Free fermionic model at the LHC, tan/3 = 40 and g z = 0.2 
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Figure 8: Free fermionic model at the LHC, tan (3 = 40 and g z = 0.1 

compared to the case of Mz> = 1.2 TeV. A similar analysis of the shapes of the resonances 
is shown in Fig. [T2] where we have chosen but this time we have varied the strength of the 
new coupling in order to show the widening of the width, which may easy the detection 
of the new neutral currents. As shown in Tab. (DO), only at large values of the couplings 
the size of the width is such to ensure a more direct identification of the resonance, which 
should probably be around 30 GeV or more, in order not to be missed. We conclude this 
section with the discussion of some results concerning the study of the variation of the 
cross section da/dQ(Q = Mz>) (on the peak) as we vary the factorization scale. In Fig. 
[13] the scale \if has been varied in the interval l/2Mz> < /if < 2Mz> for a mass Mz> = 600 
GeV. These variations are rather small over all the energy interval that we have analyzed 
and show consistently the reduction of the scale dependence of the result moving from 
LO to NLO and NNLO. The cross section is sizeable in particular above the 4 TeV scale, 
especially for larger couplings, although the presence of the resonance is not resolved in 
this figure given the small width. Finally, in Fig. [TUwe plot the total cross section as a 
function of the energy for 3 values of the new gauge couplings for Mz> = 1-2 TeV. Also in 
this case the rise of the cross section gets sizeable for larger value of the couplings. 
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Figure 9: Free fermionic model at the LHC, tan/? 
SM results through the same perturbative orders. 
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4.3 NLO/NNLO comparisons and relative differences 

We have included a set of tables which may be useful for actual experimental searches 
and comparisons. In table [5] we show the LO and in table [3] the NLO results for the 
invariant mass distributions for the first choice (800 GeV) of the mass of the extra Z' in 
all the models, and the corresponding value also for the SM. In all the cases the proximity 
among the various determinations is quite evident, except on the resonance, where the 
values show a wide variability. The pattern at NNLO, shown in table H] is similar, and the 
changes in the cross sections from NLO to NNLO in most of the cases are around 3 % or 
less. These changes are of the same order of those obtained by a study of the K-factors 
in the case of the Z resonance [20J. Also for this kinematical region, as on the Z peak 
[20] . the changes from LO to NLO are around 20-30 %, and cover the bulk of the QCD 
corrections. The last several tables describe the relative differences between the results of 
the various models and the SM, normalized to the SM values, at the various perturbative 
orders and for 3 values of the coupling constants gz = 0.05, 0.1 and 0.2. They give an 
indication of the role played by the changes in the coupling on the behaviour of these 
observables at the tails of the resonance region. In tables and [U] the region that we 
explore is between 1 and 1.5 TeV. It is rather clear from these results that for a weakly 
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NNLO FF model, M z = 1.2 TeV^ g = 0.1, tanp = 40 

NLO 
LO 

LO SM 

NLO SM 




coupled Z' (gz = 0.05) the NLO and NNLO variations respect to the SM result are 
essentially similar. The differences at NLO between the various models and the NLO SM 
are a fraction of a percent. Therefore, NNLO QCD corrections will not help in this region 
for such weakly coupled extra Z' . The differences are not more sizeable as we increase the 
new gauge coupling to 0.1, as shown in [7] and [HJ Both at NLO and NNLO the difference 
between the SM background and all the other models is smaller than 1 %. Things are not 
much better for a value of the coupling constant equal to 0.2. The differences between the 
SM and various models in this region of fast fall-off can be of the order of only 2 %, and 
just for one model ( U B — L"). Given also the small size of these cross sections, which are 
of the order of 3 x 10~ 2 fb, it is hard to separate the various contributions. Naturally, the 
situation will improve considerably if we allow a larger gauge coupling since the differences 
between signal and background can become, in principle, quite large. 



5 Conclusions 



We performed a preliminary comparative analysis of the behaviour of several models 
containing extra neutral currents in anomaly-free constructions and we discussed the 
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Figure 11: Free fermionic model and the corresponding SM results at all the three orders 
for M z > = 2.5 TeV. 



implications of the results for actual experimental searches at the LHC. Compared to 
other studies, our objective has been to compare signal and QCD background in a series 
of models, with the highest accuracy, which can be systematically performed through 
NNLO. As expected, the critical parameters in order to be able to see a signal of these 
new interactions at the new collider are the size of the gauge coupling and the mass of 
the extra gauge boson, while the specific charge assignments of the models play a minor 
role. Other parameters such as tan j3 also do not play any significant role in these types 
of searches. It is reasonable to believe that much of the potentiality for discovering the 
new resonance, if found, is its width, and all the models analyzed so far show very similar 
patterns, with a gauging of U B — L" being the one that has a slightly wider resonant 
behaviour. Being the coupling so important in order to identify which model has better 
chances to be confirmed or ruled out, it is necessary, especially in bottom- up constructions, 
to rely on more precise investigations of possible scenarios for the running of the couplings, 
which are not addressed in approaches of these types. In the case of the free fermionic U(l) 
that we have analyzed, the possibility to include these models in a more general scenario 
is natural, since they are naturally produced by a unification scheme, but is left for future 
studies. On the other hand, in these and similar models obtained either in the string 
picture or in Grand Unification, the decoupling of part of the "extra stuff" that would 
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Figure 12: Free fermionic model and the corresponding SM results at NNLO for Mz> = 2.5 
TeV for different values of g z larger than g z = 0.1. 



complicate the scenario that we have analyzed, requires extra assumptions, which would 
also affect the running of the couplings of the extra U(l)'s. These assumptions would 
introduce various alternatives on the choice of the symmetry breaking scales, threshold 
enhancements, and so on, which amount, however, to important phenomenological details 
which strongly affect this search. 

Since the V-A structure of the couplings exhibits differences with respect to other Z' 
models a measurement of forward-backward asymmetries and/or of charge asymmetries 
could be helpful [18J, but only if the gauge coupling is sizeable. The discrimination among 
the various models remains a very difficult issue for which NNLO QCD determinations, 
at least in leptoproduction, though useful, do not seem to be necessary in a first analysis. 
For those values of the mass of the extra Z' that we have considered these corrections 
cannot be isolated, while the NLO effects remain important. 
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Figure 13: Study of the fip scale dependence in the total cross section for the U(1)b-l 
model with Mz> = 0.6 TeV and g z = 0.1. Here we have chosen Mz> = Q for semplicity. 
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Figure 14: Total cross section for the Free fermionic model at NLO for three different 
values of g z and for Mz> = 1.2 TeV. Here we have chosen /ip = /ir = Q for semplicity 
and we have integrated the mass invariant distribution on the interval Mz> ± 3r^/. 
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r Mz ,(g z ) [GcV] 


9z 


M z > = 0.8 TcV 


M z > = 1.2 TeV 


M z , = 2.5 TcV 


0.02 


0.004 


0.005 


0.012 


0.05 


0.024 


0.036 


0.075 


0.1 


0.097 


0.146 


0.303 


0.2 


0.388 


0.584 


1.215 


0.3 


0.875 


1.314 


2.735 


0.4 


1.555 


2.336 


4.863 


0.5 


2.430 


3.650 


7.598 


0.6 


3.500 


5.256 


10.94 


0.7 


4.764 


7.154 


14.89 


0.8 


6.223 


9.344 


19.45 


0.9 


7.876 


11.82 


24.61 


1 


9.723 


14.60 


30.39 



Table 1: Dependence of the total width on the coupling constant g z for the free fermionic 
model with M z , = 800 GeV, M z , = 1.2 TeV and M z , = 2.5 TeV. 



25 



Acknowledgements 

We thank Simone Morelli for discussions and for various cross-checks in the numerical 
analysis. M.G. thanks the Theory Division at the University of Liverpool for hospitality 
and the Royal Society for financial support. The work of C.C. was supported (in part) 
by the European Union through the Marie Curie Research and Training Network "Uni- 
versenet" (MRTN-CT-2006-035863) and by The Interreg II Crete-Cyprus Program. He 
thanks the Theory group at Crete for hospitality. The work of A.E.F. is supported in 
part the STFC. 



26 



da LO /dQ [pb/GcV], M' z = 800, g z = 0.1, tan/3 = 40, Candia cvol. 


Q [GcV] 


(Jlo(Q) FFM 


<?lo{Q) U(1) B -L 


vlo{Q) U(l) q+U 


vlo(Q) f/(l)io+5 


<Jlo(Q) U(l) d - U 


<Jlo(SM) 


750 


1.1101 • 10~ 4 


1.0854- 10~ 4 


1.0854- 10" 4 


1.1017- 10~ 4 


1.1011 • 10~ 4 


1.1033- 10~ 4 


761 


1.0355- 10- 4 


1.0050- 10~ 4 


1.0050- 10- 4 


1.0250- 10~ 4 


1.0243 • 10~ 4 


1.0269- 10~ 4 


773 


9.6852 • 10- 5 


9.2759 • 10~ 5 


9.2759 • 10~ 5 


9.5421 • 10" 5 


9.5315 • 10~ 5 


9.5674 • 10~ 5 


784 


9.1225- 10- 5 


8.4635 • 10~ 5 


8.4635 • 10~ 5 


8.8822 • 10" 5 


8.8633 • 10~ 5 


8.9212 • 10~ 5 


796 


9.1654- 10" 5 


7.2110- 10~ 5 


7.2111 • 10~ 5 


8.2428- 10" 5 


8.1409 • 10~ 5 


8.3259 • 10~ 5 


800 


1.6448- 10~ 2 


4.2388- 10~ 2 


2.3928- 10" 2 


1.9570- 10~ 2 


4.3085- 10- 2 


8.1086- 10- 5 


800 


4.9572 • 10~ 4 


1.9334- 10~ 3 


1.8812-10- 3 


2.8631 • 10~ 4 


1.5888- 10- 4 


8.0955 • 10~ 5 


801 


1.7452- 10~ 4 


6.8269 • 10~ 4 


6.7771 • 10~ 4 


1.4480- 10~ 4 


1.1071 • 10~ 4 


8.0839 • 10~ 5 


839 


6.4010- 10" 5 


6.6355 • 10~ 5 


6.6355 • 10~ 5 


6.4761 • 10" 5 


6.4800 • 10~ 5 


6.4607 • 10~ 5 


868 


5.4301 • 10" 5 


5.5480 • 10~ 5 


5.5480 • 10~ 5 


5.4686 • 10" 5 


5.4706 • 10~ 5 


5.4610 • 10~ 5 


900 


4.5656- 10" 5 


4.6371 • 10" 5 


4.6371 • 10~ 5 


4.5892 • 10" 5 


4.5904 • 10~ 5 


4.5847 • 10~ 5 



Table 2: LO invariant mass distributions 



da NLO /dQ [pb/GcV], M' z = 800, g z = 0.1, tan/3 = 40, Candia evol. 


Q [GcV] 


>Jnlo(Q) FFM 


onlo{Q) U(1) B _ L 


<Jnlo(Q) U(l) q+U 


ctjvlo(Q) t/(l)io+5 


vnlo(Q) U(l) d _ u 


(Jnlo{SM) 


750 


1.4362- 10~ 4 


1.4048- 10~ 4 


1.4048 • 10~ 4 


1.4257- IO" 4 


1.4248- 10~ 4 


1.4276- 10~ 4 


761 


1.3394- 10~ 4 


1.3008- IO" 4 


1.3008 • 10~ 4 


1.3263- IO" 4 


1.3252- 10~ 4 


1.3287- 10~ 4 


773 


1.2526- 10~ 4 


1.2006- 10~ 4 


1.2006- IO" 4 


1.2346 • 10~ 4 


1.2331 • 10~ 4 


1.2377- 10~ 4 


784 


1.1794- 10~ 4 


1.0958- IO" 4 


1.0958- IO" 4 


1.1492- IO" 4 


1.1465- 10~ 4 


1.1540- 10~ 4 


796 


1.1834- 10~ 4 


9.3647- 10" 5 


9.3648 • 10~ 5 


1.0671 • 10~ 4 


1.0530- IO" 4 


1.0769- IO" 4 


800 


2.1411 • 10~ 2 


5.4992 • 10~ 2 


3.1043-10" 2 


2.5389- 10~ 2 


5.5896- IO" 2 


1.0487- 10~ 4 


800 


6.4893- IO" 4 


2.5014- 10~ 3 


2.4339 • IO" 3 


3.6947- IO" 4 


2.0566- 10~ 4 


1.0470- 10~ 4 


801 


2.2882- 10~ 4 


8.8180- IO" 4 


8.7537- IO" 4 


1.8666- IO" 4 


1.4323- 10~ 4 


1.0455- 10~ 4 


839 


8.2772 • 10~ 5 


8.5749 • 10~ 5 


8.5749 • 10~ 5 


8.3712 • 10~ 5 


8.3771 • 10~ 5 


8.3523- IO" 5 


868 


7.0183- IO" 5 


7.1679- IO" 5 


7.1679 • 10~ 5 


7.0666 • 10~ 5 


7.0696 • 10~ 5 


7.0573 • 10~ 5 


900 


5.8982 • 10~ 5 


5.9888 • 10~ 5 


5.9888 • 10~ 5 


5.9278 • 10~ 5 


5.9296 • 10~ 5 


5.9222 • 10~ 5 



Table 3: NLO distributions for 750 < Q < 900 GeV 



Acjnnlo/AQ [pb/GeV], M' z = 800, g z = 0.1, tan/3 = 40, Candia evol. 


Q [GcV] 


(Jnnlo(Q) FFM 


& N N Lo{Q) U(1)b-L 


O'NNLoiQ) U(l) q+U 


&NNLo(Q) C^(l)lO+5 


cnnlo(Q) U(l)d- U 


(Jnnlo(SM) 


750 


1.4793- 10~ 4 


1.4472 • 10~ 4 


1.4472 • 10~ 4 


1.4686- 10~ 4 


1.4676 • 10~ 4 


1.4705 • 10~ 4 


761 


1.3803- 10^ 4 


1.3407- IO" 4 


1.3407- 10~ 4 


1.3669- IO" 4 


1.3657- IO" 4 


1.3693- IO" 4 


773 


1.2914- 10~ 4 


1.2382 • 10~ 4 


1.2382- 10~ 4 


1.2730- 10~ 4 


1.2714- IO" 4 


1.2762- 10~ 4 


784 


1.2164- IO" 4 


1.1308- IO" 4 


1.1308- 10~ 4 


1.1856- IO" 4 


1.1827- IO" 4 


1.1904- 10~ 4 


796 


1.2207- IO" 4 


9.6772 • 10~ 5 


9.6773 • 10~ 5 


1.1017- 10~ 4 


1.0867- IO" 4 


1.1114- 10~ 4 


800 


2.2140- 10~ 2 


5.6805 • 10~ 2 


3.2066- IO" 2 


2.6233- 10~ 2 


5.7755- IO" 2 


1.0825- 10~ 4 


800 


6.7222- 10~ 4 


2.5818-10" 3 


2.5121 • IO" 3 


3.8114- IO" 4 


2.1235 • 10~ 4 


1.0808- 10~ 4 


801 


2.3717- 10~ 4 


9.0971 • 10~ 4 


9.0307- IO" 4 


1.9249- 10~ 4 


1.4787- IO" 4 


1.0793- IO" 4 


839 


8.5581 • 10~ 5 


8.8638 • 10~ 5 


8.8638 • 10~ 5 


8.6542 • 10~ 5 


8.6606 • 10" 5 


8.6349 • 10" 5 


868 


7.2645 • 10~ 5 


7.4182 • 10~ 5 


7.4182 • 10~ 5 


7.3139- IO" 5 


7.3171 • 10" 5 


7.3044 • 10" 5 


900 


6.1122- 10~ 5 


6.2054 • 10~ 5 


6.2054- 10~ 5 


6.1425- 10~ 5 


6.1444- 10" 5 


6.1368 • 10" 5 



Table 4: NNLO distributions for 750 < Q < 900 GeV 



\°n™ ~ <io\lolio% i M 'z = 80 °: 9z = 0.05, tan = 40, Candia evol. 


Q [GeV] 


^(Q)[pb/GeV] 


A FFMo/ 










1000 


3.5146 • 10" 5 


6.5325- 10~ 2 


1.6003- lO" 1 


1.6003- 10" 1 


1.0162- 10" 2 


1.4126- 10" 2 


1015 


3.2618 • 10" 5 


6.2528- 10" 2 


1.5220- 10" 1 


1.5220- 10" 1 


9.5155- 10" 3 


1.3203- 10" 2 


1030 


3.0299 • 10" 5 


6.0105 • 10~ 2 


1.4541 • 10" 1 


1.4541 • 10" 1 


8.9574-10" 3 


1.2402 • 10" 2 


1045 


2.8168 • 10" 5 


5.7987- 10" 2 


1.3947- 10" 1 


1.3947- 10" 1 


8.4716 • 10~ 3 


1.1701 • 10~ 2 


1060 


2.6209 • 10" 5 


5.6121 • 10~ 2 


1.3424- 10" 1 


1.3424- 10" 1 


8.0455 • 10~ 3 


1.1083 • 10" 2 


1165 


1.6156- 10" 5 


4.7509- 10- 2 


1.1003- 10- 1 


1.1003- 10" 1 


6.1091 • 10^3 


8.2061 • 10~ 3 


1210 


1.3265 • 10" 5 


4.5240- 10" 2 


1.0361 • 10" 1 


1.0361 • 10" 1 


5.6126 • 10~ 3 


7.4369 • lO- 3 


1250 


1.1183- 10" 5 


4.3636- 10- 2 


9.9050- 10- 2 


9.9050- 10- 2 


5.2687-10- 3 


6.8881 • lO" 3 


1355 


7.2763 • 10" 6 


4.0639- 10- 2 


9.0453- 10- 2 


9.0453- 10- 2 


4.6479 • 10~ 3 


5.8404- lO" 3 


1425 


5.5361 • 10~ 6 


3.9279- 10- 2 


8.6492 • 10~ 2 


8.6492- 10~ 2 


4.3830 • 10~ 3 


5.3492 • lO" 3 


1500 


4.1734- 10" 6 


3.8186- 10- 2 


8.3253- 10- 2 


8.3253- 10- 2 


4.1841 • 10^3 


4.9401 • lO" 3 



Table 5: Percentage differences at NLO. We define A^ /o = \<r%% - cr^/af^ where i = FFM, B - L, q + u, 10 + 5, d - u. 



CO 



Knio ~ <mo\/° S nnio% , M' z = 800, g z = 0.05, tan P = 40, Candia evol. 


Q [GcV] 


« (Q)[pb/GcV] 


A FFMo/ 






a10+5ck 
A nnio /o 




1000 


T 3.6546 • 10~ 5 


6.4565 • 10" 2 


1.5879- 10" 1 


1.5879- 10" 1 


9.8298- lO" 3 


1.4114- lO" 2 


1015 


3.3935 • lO" 5 


6.1789- lO" 2 


1.5099- lO" 1 


1.5099- lO" 1 


9.1914- lO" 3 


1.3191 • lO" 2 


1030 


3.1537- lO" 5 


5.9383- lO" 2 


1.4423- lO" 1 


1.4423- 10" 1 


8.6399-10-3 


1.2390- lO" 2 


1045 


2.9334 • 10~ 5 


5.7279- lO" 2 


1.3831 • lO" 1 


1.3831 • lO" 1 


8.1595-10-3 


1.1689- lO" 2 


1060 


2.7306 • 10~ 5 


5.5425 • 10~ 2 


1.3310- lO" 1 


1.3310- 10" 1 


7.7381 • lO" 3 


1.1071 • lO" 2 


1165 


1.6888- lO" 5 


4.6857- lO" 2 


1.0895- lO" 1 


1.0895- lO" 1 


5.8167-10-3 


8.1924-10-3 


1210 


1.3884- 10" 5 


4.4594- 10~ 2 


1.0254- 10" 1 


1.0254- lO" 1 


5.3213- lO" 3 


7.4223 • lO" 3 


1250 


1.1718- 10~ 5 


4.2991 • 10~ 2 


9.7986- 10" 2 


9.7986- 10" 2 


4.9768-10-3 


6.8729 • lO" 3 


1355 


7.6472 • 10" 6 


3.9988- lO" 2 


8.9375 • 10~ 2 


8.9375- lO" 2 


4.3502- lO" 3 


5.8236-10-3 


1425 


5.8293- lO" 6 


3.8618- lO" 2 


8.5395- lO" 2 


8.5395- lO" 2 


4.0792 • lO" 3 


5.3312-10-3 


1500 


4.4031 • 10~ 6 


3.7510- lO" 2 


8.2129- 10~ 2 


8.2129- 10~ 2 


3.8720-10-3 


4.9211 • lO" 3 



Table 6: Percentage differences at NNLO. We define A l nnlo = \cr^ — cr l nnlo \/cr. 



- <lo\/°™% . M z = 800 > 5, = 0.1, tan/3 - 40, Candia evol. 




„SM ( r\ \ r^u /ri^Mi 
a nlo (Q)lP b /GcVJ 


A FFM 0/ 


A nlo % 


A nlo % 


A 10+5o/ 


A d — uoy 


1000 


3.5146 • 10" 5 


2.4555 • 10" 1 


6.5950 • 10 _1 


6.5950 • 10 _1 


5.5268 • 10~ 2 


7.0677 • 10~ 2 


1015 


3.2618 • 10~ 5 


2.3454 • 10 _1 


6.2764 • 10 _1 


6.2764 • 10 _1 


5.2612 • 10~ 2 


6.6957- 10~ 2 


1030 


3.0299 • 10~ 5 


2.2500 • 10" 1 


6.0005 • 10 _1 


6.0005 • 10 _1 


5.0320 • 10~ 2 


6.3730 • 10~ 2 


1045 


2.8168 • 10~ 5 


2.1666 • 10 _1 


5.7593- 10" 1 


5.7593- 10" 1 


4.8324 • 10~ 2 


6.0904- 10~ 2 


1060 


2.6209 • 10" 5 


2.0931 • 10" 1 


5.5469 • 10 _1 


5.5469 • 10 _1 


4.6574 • 10~ 2 


5.8410 • 10" 2 


1165 


1.6156 • 10~ 5 


1.7534- 10" 1 


4.5649 • 10 _1 


4.5649 • 10 _1 


3.8610 • 10~ 2 


4.6791 • 10~ 2 


1210 


1.3265 • 10~ 5 


1.6639 • 10 _1 


4.3048 • 10 _1 


4.3048 • 10 _1 


3.6562 • 10~ 2 


4.3677 • 10~ 2 


1250 


1.1183- 10~ 5 


1.6007- lO" 1 


4.1202 • 10" 1 


4.1202- 10" 1 


3.5137- 10" 2 


4.1449 • 10~ 2 


1355 


7.2763 • 10~ 6 


1.4826- 10" 1 


3.7722- 10" 1 


3.7722- lO" 1 


3.2554- 10" 2 


3.7189- 10- 2 


1425 


5.5361 • 10~ 6 


1.4291 • lO" 1 


3.6118- 10" 1 


3.6118- 10" 1 


3.1440 • 10~ 2 


3.5182 • 10~ 2 


1500 


4.1734- 10~ 6 


1.3861 • 10- 1 


3.4806- 10- 1 


3.4806- 10- 1 


3.0591 • 10~ 2 


3.3504- 10~ 2 



Table 7: Percentage differences at NLO for g z = 0.1. Here and in the following we use the same notation of the previous tables. 



CO 

o 



Wnnio ~ <nio\/Vnnio% > M 'z = 800 > 9z = 0.1, tan/3 = 40, Candia evol. 


Q [GcV] 


<^ (Q)[pb/GeV] 


A FFM 01 






Al0+5 % 




1000 


3.6546 • 10~ 5 


2.4248 • lO" 1 


6.5458- 10- 1 


6.5458- 10- 1 


5.3970- 10- 2 


7.0662 • 10~ 2 


1015 


3.3935 • 10~ 5 


2.3155- 10- 1 


6.2284- lO" 1 


6.2284- 10" 1 


5.1346- 10~ 2 


6.6942 • 10~ 2 


1030 


3.1537- 10- 5 


2.2208- lO" 1 


5.9536- 10- 1 


5.9536- 10- 1 


4.9081 • 10~ 2 


6.3715- 10- 2 


1045 


2.9334- 10~ 5 


2.1379- 10" 1 


5.7133- 10" 1 


5.7133- 10" 1 


4.7108- 10~ 2 


6.0888- 10" 2 


1060 


2.7306 • 10~ 5 


2.0649- 10" 1 


5.5017- 10- 1 


5.5017- 10- 1 


4.5376- 10- 2 


5.8394- 10- 2 


1165 


1.6888 • 10~ 5 


1.7269- lO" 1 


4.5223- lO" 1 


4.5223- lO" 1 


3.7478- 10~ 2 


4.6774- 10" 2 


1210 


1.3884- 10~ 5 


1.6377- 10- 1 


4.2627- 10" 1 


4.2627- 10" 1 


3.5438- 10- 2 


4.3659- 10- 2 


1250 


1.1718- 10~ 5 


1.5744- 10" 1 


4.0781 • lO" 1 


4.0781 • lO" 1 


3.4013- 10~ 2 


4.1431 • 10~ 2 


1355 


7.6472 • 10~ 6 


1.4560- lO" 1 


3.7296- 10- 1 


3.7296- 10- 1 


3.1411 • 10~ 2 


3.7169- 10- 2 


1425 


5.8293- 10~ 6 


1.4021 • 10" 1 


3.5685- 10- 1 


3.5685- 10- 1 


3.0275- 10- 2 


3.5160- 10- 2 


1500 


4.4031 • 10~ 6 


1.3585- 10- 1 


3.4362- 10" 1 


3.4362- 10" 1 


2.9397- 10~ 2 


3.3481 • 10~ 2 



Table 8: Percentage differences at NNLO for g z — 0.1 



knto - <lo\/°™% . M z = 800 > fl* = °- 2 > tan/ 3 = 40 > Candia evoL 






A FFMW 


A B — Lq/ 


A <?+«(>/ 


A 10+5o/ 
A nlo % 


A d — UQ/ 


1000 


3.5146 • 10~ 5 


9.4061 • lO^ 1 


2.7377 • 10+° 


2.7377 • 10+° 


2.4462 • 10" 1 


2.9911 • 10" 1 


1015 


3.2618 • 10~ 5 


8.9927- 10" 1 


2.6020 • 10+° 


2.6020 • 10+° 


2.3306- 10" 1 


2.8399- 10" 1 


1030 


3.0299 • 10~ 5 


8.6334 • 10 _1 


2.4847 • 10+° 


2.4847- 10+° 


2.2312 • 10" 1 


2.7088 • 10 _1 


1045 


2.8168 • 10~ 5 


8.3183 • 10" 1 


2.3825 • 10+° 


2.3825 • 10+° 


2.1448 • 10 _1 


2.5940 • 10 _1 


1060 


2.6209 • 10~ 5 


8.0401 • 10" 1 


2.2926 • 10+° 


2.2926 • 10+° 


2.0693 • 10 _1 


2.4928 • 10" 1 


1165 


1.6156 • 10~ 5 


6.7477- 10" 1 


1.8795 • 10+° 


1.8795 • 10+° 


1.7274 • 10 _1 


2.0217- 10 _1 


1210 


1.3265 • 10~ 5 


6.4052 • 10 _1 


1.7707 • 10 +0 


1.7707 • 10 +0 


1.6400 • 10 _1 


1.8955 • 10 _1 


1250 


1.1183- 10~ 5 


6.1627- 10" 1 


1.6937- 10+° 


1.6937- 10+° 


1.5792- 10" 1 


1.8052- 10" 1 


1355 


7.2763- 10~ 6 


5.7094- 10- 1 


1.5487- 10+° 


1.5487- 10+° 


1.4689- 10" 1 


1.6326- 10" 1 


1425 


5.5361 • 10~ 6 


5.5039 • 10" 1 


1.4820 • 10+° 


1.4820- 10+° 


1.4212 • lO" 1 


1.5512 • 10" 1 


1500 


4.1734- 10~ 6 


5.3389- 10- 1 


1.4275 • 10+° 


1.4275 • 10+° 


1.3846- 10- 1 


1.4832- 10" 1 



Table 9: Percentage differences at NLO for g z = 0.2 



Wnnio ~ <nio\/Vnnio% > M 'z = 800 > 9* = °- 2 , tan I 3 = 40 > Candia evol. 


Q [GcV] 


<7^ (Q)[pb/GeV] 


A FFMov 


A Z^% 


A l + n l% 


A 10+5 % 
^nnlo /0 




1000 


3.6546 • 10~ 5 


9.2821 • 10" 1 


2.7182- 10+° 


2.7182- 10+° 


2.3947- 10" 1 


2.9909- 10" 1 


1015 


3.3935 • 10~ 5 


8.8720- 10- 1 


2.5829 • 10+° 


2.5829- 10+° 


2.2804- lO" 1 


2.8397- 10- 1 


1030 


3.1537- 10- 5 


8.5154- lO" 1 


2.4661 • 10+° 


2.4661 • 10+° 


2.1821 • 10" 1 


2.7085- 10- 1 


1045 


2.9334 • 10~ 5 


8.2026- 10" 1 


2.3642 • 10+° 


2.3642 • 10+° 


2.0966- lO" 1 


2.5938- 10" 1 


1060 


2.7306 • 10~ 5 


7.9263- 10- 1 


2.2747- 10+° 


2.2747- 10+° 


2.0218- 10" 1 


2.4925- 10" 1 


1165 


1.6888- 10" 5 


6.6409- 10- 1 


1.8626 • 10+° 


1.8626- 10+° 


1.6826- lO" 1 


2.0214- 10" 1 


1210 


1.3884- 10~ 5 


6.2993- 10- 1 


1.7540 • 10+° 


1.7540- 10+° 


1.5955- 10- 1 


1.8952- 10- 1 


1250 


1.1718- 10~ 5 


6.0571 • 10- 1 


1.6769 • 10+° 


1.6769- 10+° 


1.5347- 10- 1 


1.8049- 10" 1 


1355 


7.6472 • 10~ 6 


5.6024- lO" 1 


1.5317- 10+° 


1.5317- 10+° 


1.4237- 10" 1 


1.6323- lO" 1 


1425 


5.8293- 10~ 6 


5.3951 • 10- 1 


1.4647- 10+° 


1.4647- 10+° 


1.3751 • 10- 1 


1.5509- 10- 1 


1500 


4.4031 • 10~ 6 


5.2277- 10- 1 


1.4098 • 10+° 


1.4098- 10+° 


1.3374- 10" 1 


1.4828- 10" 1 



Table 10: Percentage differences at NNLO for g z = 0.2 



da nnlo /dQ [pb/GeV] for the FF model with M z , = 2.5 TeV, tan/3 = 40, Candia evol. 



Q [TcV] 


9z = 0.1 


t _ 5,-0.4 


9z = 0.6 


g z = 0.8 


9z = 1 




2.400 


2.6475 • 10~ 7 


3.3941 • 10~ 7 


5.0947- 10- 7 


8.7995 • 10~ 7 


1.5720- 10~ 6 


2.6141 • 10~ 7 


2.423 


2.4961 • 10~ 7 


3.5212 • 10~ 7 


6.0291 • 10~ 7 


1.1654- 10" 6 


2.2223- 10~ 6 


2.4543 • 10~ 7 


2.446 


2.3629- 10~ 7 


4.0068 • 10~ 7 


8.4077- 10- 7 


1.8529- 10" 6 


3.7317- 10- 6 


2.3050- 10- 7 


2.469 


2.2656- 10~ 7 


6.0047- 10- 7 


1.7162 • 10~ 6 


4.2536- 10" 6 


8.5322- 10~ 6 


2.1654- 10~ 7 


2.492 


2.4932- 10~ 7 


3.7446 • 10~ 6 


1.2697- 10" 5 


2.3281 • 10~ 5 


3.0409 • 10~ 5 


2.0349 • 10~ 7 


2.5000 


3.7947 • 10~ 5 


3.7947 • 10~ 5 


3.7947 • 10~ 5 


3.7947 • 10~ 5 


3.7947 • 10~ 5 


1.9900- 10- 7 


2.5003 


8.5271 • 10~ 6 


3.7283 • 10~ 5 


3.7757 • 10~ 5 


3.7858 • 10~ 5 


3.7892 • 10~ 5 


1.9886- 10- 7 


2.5005 


2.7949 • 10~ 6 


3.5983 • 10~ 5 


3.7438 • 10~ 5 


3.7730 • 10~ 5 


3.7824 • 10~ 5 


1.9873- 10" 7 


2.5770 


1.5907- 10- 7 


1.4769 • 10~ 7 


2.2368- 10" 7 


5.0120- 10- 7 


1.1340- 10~ 6 


1.6192- 10~ 7 


2.636 


1.3692- 10- 7 


1.2412 • 10~ 7 


1.3364- 10" 7 


1.9772 • 10~ 7 


3.6561 • 10~ 7 


1.3839- 10- 7 


2.700 


1.1628- 10~ 7 


1.0680- 10" 7 


1.0536- 10~ 7 


1.2481 • 10~ 7 


1.8637- 10" 7 


1.1718- 10" 7 



Table 11: NNLO cross sections for the FF model with a Mz> = 2.5 TeV for values of the 
coupling constant g z larger than g z = 0.1 
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